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ABSTRACT 

A detailed analysis of NGST multiband photometry applied to the reference case 
of the study of high-redshift galaxies has been carried out with simulations based on 
galaxy SEDs derived from the currently available empirical and model templates and 
on plausible standard filter-sets. 

In order to correctly identify star forming galaxies in the redshift range 5 < z < 9 and 
early-type galaxies above z > 2 and avoid confusion with other SEDs it is mandatory 
to have photometric information in optical bands, besides a standard IR filter-set like 
FllO F160 K L M. In particular by adding the ^606, 1814 and z — Gunn filters a 
good discrimination is obtained above z > 5 for star forming galaxies and z > \ for 
early-types. The case for an extension of the NGST wavelength domain to the optical 
range is therefore strongly supported by this analysis. 

The effects of leaks in the filter blocking have also been investigated. In spite of 
rather pessimistic assumptions (a constant leak at a level of 10~^ of the peak trans- 
mission over the whole spectral range or a leak of Gaussian shape placed at 1.5 
times the effective wavelength of the filter with an amplitude of 10~^ of the peak 
transmission and a width of 5% of the central wavelength) the effects are not dra- 
matic: the accuracy of the determination of photometric redshifts with a standard 
y606 7814 z — Gunn FllO F160 K L M filter set does not significantly deteriorate for 
a sample limited to Mab ~ 30. In the range 5.5 < z < 12.5 a rather high accuracy with 
a typical error \zsim — Zphot\ < 0.2 is achieved. 
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With the filter blocking achieved by the present standard technology, extended and 
continuous spectral coverage appears to be the driving factor to maximize the scientific 
output. 

1. Introduction 

The NGST Design Reference Mission (DRM) is a representative science program elaborated in 
sufficient detail to aid in the development of functional requirements for the NGST mission. One of 
the main results of the DRM has been the identification of a large fraction of subprograms (8 out 
of 12) requiring observations at optical wavelengths. Such an extension poses significant technical 
challenges. In particular it requires an efficient blocking of photons at undesired wavelengths over 
a range covering more than three octaves. In the following we examine the consequences of leaks 
in the transmission of filters for a typical programme of determination of photometric redshifts of 
high-redshift galaxies. 

2. Baseline of the Instrument 

Table 1 shows the main characteristics of the Visible-NIR camera as conceived by the Visible 
and NIR camera Subcommittee of the "Ad Hoc Science Working Group" (ASWG). Two options 
are envisaged: an optimal and a bare bones version. More details can be found at the URL 
http : //www . stecf . org/ngst/ instruments . html. 

3. Band-passes of the filters 

We have assumed as reference bandpasses for the filters of the optical-NIR camera the standard 
V606 1814 z - Gunn FllO F160 K L M set commonly used on HST with WFPC2 (1^606 7814) or 
NICMOS (FllO F160) and from the ground (Fig. 1). Details about the filter properties are given 



Table 1. Visible-NIR Camera 





Optimal 


Bare Bones 


A range 


0.6 — 5//m 


1.0 - 5.0/xm 


Sampling 


0.03" 


0.03" 


FOV 


4' X 4' 


2' X 2' 


A/AA 


< 10 


< 10 
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in Table 2 

The actual choice of the filters to be used in space with NGST may be different but the 
arguments developed in the following do not depend critically on this detail. 

4. Typical SEDs of the sources and the locus of galaxies in multicolor space 

To simulate a range of galaxy SEDs we have adopted the templates of Coleman et al. (1980) 

for a typical elliptical, Sbc, Scd and irregular galaxy plus a GISSEL (Bruzual and Chariot 2001) 
template for a very blue galaxy (with a constant star formation rate and an age of 0.1 Gyr). The 
SEDs are shown in Fig. 2 and in the following will be referred to as extended CWW templates. 

Fig. 3 shows the evolution of the apparent magnitude of typical galaxies of IO^^Mq in the 
-F606 and L bands and gives a rough idea of the magnitudes of interest. It can be appreciated 
how an elliptical fades away relatively quickly in the F606 band, while other star-forming galaxies 
remain - as observed - significantly brighter. 

The upper panel of Fig. 4 shows the evolution of the colors of the galaxies described in the 
previous section in a (F606 — F814) vs. (FllO — K) plane. The color tracks are limited to z < 3.5 
for an elliptical SED and to z < 5 for Sbc galaxies. With this combination of bandpasses it is 
straightforward to identify star- forming galaxies with z > 5 and ellipticals at z > 1. The lower 
panel shows the same in a (FllO — F160) vs. (F160 — K) plane and demonstrates how, without 
optical bands, it becomes difficult to isolate star-forming galaxies at z < 9. 

5. Photometric redshifts and the effects of "(red) leaks" 

In order to quantify the different performance of different sets of filters and the effects of 
"red-leaks" in a filter configuration extending from the visible to the infrared we have carried out 
simulations of the assignment of photometric redshifts to a sample of faint galaxies. To this end, 
we have tried to generate a simulated photometric catalog covering a very large range of galaxy 
spectral properties and to estimate photometric redshifts from an economic and independent a set 
of templates as possible. Galaxies were simulated in a redshift interval < z < 20 according to 
GISSEL library. This spectral synthesis model is controlled by a number of free parameters. The 
star formation rate for a galaxy at a given age is governed by the assumed e — folding star formation 
time-scale r. Several values of r and galaxy ages are necessary to reproduce the different observed 
spectral types. For the shape of the initial mass function (IMF) we have assumed a Salpeter law, 
which turns out not to be a critical choice. In addition to the GISSEL parameterization, we have 
added the internal reddening of each galaxy by applying the observed attenuation law for local 
starburst galaxies by Calzetti et al. (1997), with the Eb-v normalization varying between and 
0.2 mag. We have also included the Lyman absorption produced by the intergalactic medium as 
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Table 2. Properties of the assumed filter set 



Name 


< A > 


FWHM 


AB -Yega 




(A) 


(A) 


correction 


F606 


6031 


2050 


0.096 


F814 


8011 


1400 


0.417 


z-Gunn 


9125 


1200 


0.527 


FllO 


11288 


5239 


0.698 


F160 


16096 


3750 


1.313 


Ks 


21640 


2750 


1.841 


L 


37789 


5800 


2.917 


M 


47657 


6400 


3.380 



F814 F110 Ks M 



0.8 



6 



F60 

■ i\ 


ez^unn 

LA 










'L 1 








11 


II J 


I 


u 




i J 


1 









lO-i 2x10-* 3x10-* 4x10" 5x10" 



A(A) 

Fig. 1. — The reference filter set used for the simulations. 
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Fig. 2. — The spectral energy distribution of five template galaxies. 
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Fig. 4. — Redshift evolution of the colors of galaxies. Upper panel: (-F606 — -F818) vs. (FllO — K). 
Lower panel: (FllO - F160) vs. (F160 - K). 
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a function of redshift. More details on this type of simulations can be found in Arnouts et al. 
(1999). Fig. 5 shows the galaxy counts as a function of the apparent magnitude, the age and 
redshift distribution, the z — mag plane characterizing the simulation. 

Magnitudes and colors have been derived from the simulated SEDs by convolution with the 
filter passbands described in Sect. 3. Measurement errors have been added, according to the NGST 
mission simulator (http://www.ngst.stsci.edu/nms/niain/), 
assuming an exposure time of lO^s. 

From this set of "simulated observed colors" photometric redshifts have been derived using a 
standard procedure of optimization, comparing the observed fluxes with errors to the extended 
CWW set of templates described in Sect. 4. 

The has been computed as: 



where Fobserved,i is the flux observed in a given filter i, Oi is its uncertainty, Ftemplate,i is the flux 
of the template in the same filter and the sum is over the various filters. The template fiuxes have 
been normalized to the observed ones with the factor A minimizing the value : 



Fig. 6 shows the comparison between photometric redshift estimates obtained with a V^GOG /814 
z — Gunn FllO F160 K L M filter set and a filter set limited to the infrared bands. As suggested 
in the previous section in order to correctly identify star forming galaxies in the redshift range 
5 < z < 9 and avoid confusion with other SEDs it appears mandatory to have, besides the 
IR filters, photometric information in the optical bands. In the following we will adopt the set 
F606 7814 z - Gunn FllO F160 K L M as standard. 

In order to explore the effects of filter leaks the following approach was adopted: 

1. two types of leaks have been foreseen: a constant leak at a level of 10^^ of the peak trans- 
mission over the whole spectral range and a leak of Gaussian shape placed at 1.5 times the 
effective wavelength of the filter with an amplitude of 10~^ of the peak transmission and a 
width of 5% of the central wavelength. The effect of Gaussian leaks on the filter transmission 
is shown in Fig. 9 

2. The leaks have been included in the system spectral response and new magnitudes and colors 
have been produced with the GISSEL models. 

3. The new colors have been used to produce photometric redshifts with the extended CWW 
templates. No knowledge of the leaks was used in this step, i.e. the system responses assumed 
in the z-phot estimation are those of Fig. 1, without leaks. 




(1) 




F'2 



(2) 
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Fig. 5. — Characteristic distributions of the simulation of galaxies derived from GISSEL models: 
differential counts as a function of the apparent magnitude (upper left), redshift distribution (upper 
and lower right), the z — mag plane (lower left). In the upper left panel the continuous and dashed 
lines refer to the differential distribution in M and / magnitude, respectively. In the upper right 
panel the dashed line refers to galaxies with ages larger than 1 Gyr, the continuous line to ages 
between 10^ and 10^ yr. The colors of the dots in the lower left panel (in the electronic version) 
correspond to different galaxy types with the same color coding of Fig. 4. 
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Fig. 6. — Comparison between the accuracy of the photometric redshift estimates with a filter set 
including IR and optical bands (upper panel V606 /814 z — Gunn F110( J) FlQO{H) K L M) and 
a filter set limited to infrared bands (lower panel). The dot-dashed, long-dashed, short-dashed and 
continuous lines correspond to the 95, 90, 68 and 50 percentile of the absolute deviation \zsim— Zphot\-, 
respectively. The galaxy sample is limited to Mab < 29. 
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Fig. 7. — Typical errors of the photometric redshifts with a 1^606 1814 z — 
Gunn FllO(J) -F160(-ff) K L M filter set. The upper panel shows the result obtained assum- 
ing only the standard photometric error, the lower panel includes the effects of a constant leak at a 
level of 10~^ of the peak transmission over the whole spectral range. The dot-dashed, long-dashed, 
short-dashed and continuous lines correspond to the 95, 90, 68 and 50 percentile of the absolute 
deviation \zsim — Zphot\, respectively. The galaxy sample is limited to Mab < 28. 
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Fig. 8. — Same as Fig. 7, but down to Mab = 30. 
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Fig. 9.— The response of the 1^606 1814 z - Gunn FllO F160 K L M filter with the addition of a 
leak of Gaussian shape placed at 1.5 times the effective wavelength of the filter with an amplitude 
of 10"'^ of the peak transmission and a width of 5% of the central wavelength. 
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Fig. 10. — Typical errors of the photometric redshifts with a ^606 1814 z— Gunn FIW F160 K L M 
filter set. The upper panel shows the result obtained assuming only the standard photometric error, 
the lower panel includes the effects of a leak of Gaussian shape placed at 1.5 times the effective 
wavelength of the filter with an amplitude of 10~^ of the peak transmission and a width of 5% of the 
central wavelength. The dot-dashed, long-dashed, short-dashed and continuous lines correspond to 
the 95, 90, 68 and 50 percentile of the absolute deviation \zsim — Zphot\, respectively. The galaxy 
sample is limited to Mab < 28. 
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Fig. 11.— Same as Fig. 10 but down to Mab = 30. 
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The typical errors deriving from photometric uncertainties and unaccounted leaks are illus- 
trated in Fig. 7-8 and 10-11. The dot-dashed, long-dashed, short-dashed and continuous lines 
correspond to the 95, 90, 68 and 50 percentile of the absolute deviation \zsim — Zphot\, respectively. 
It can be seen that errors are in all cases relatively large at z < 5. This is an obvious effect due to 
the lack of bandpasses bluer than Vqqq, and the consequent difficulty to identify the Lyman break 
at relatively low redshift and distinguish it, for example, from the 4000 A break at lower redshift. 

The effect of the leaks, despite the rather pessimistic assumptions, is in general a moderate 
increase of the error level. For example for 2; = 11 and Mab < 30 the 68 percentile corresponds to 
\zsim — Zphot\ =0.23 with the standard photometric errors, while it increases at \zsim — Zphotl =0.25 
with the unaccounted constant leak at a level of 10~^ of the peak transmission over the whole 
spectral range. No qualitative change of the accuracy is observed when the leaks are introduced. 
No relevant difference is also found between the continuous and the Gaussian leak cases. Again 
for z = 11 and Mab < 30 the 68 percentile corresponds to \zsim — Zphot] = 0.27 when the leak of 
Gaussian shape is added. In the range 5.5 < z < 12.5 a rather good accuracy is achieved with a 
typical error \zsim — Zpjiot\ ^ 0.2 in all cases down to Mab ~ 30. The situation deteriorates slightly 
around z ~ 15, but the result is due mainly to the insufficient sampling of the SEDs between the 
Kg and the M filters, in particular to the "hole" at 2.9 /Ltm. 

The importance of the leaks, relative to the uncertainty due to photometric errors, decreases 
with decreasing flux. At Mab < 30 and 2; = 11 \zsim ~ Zphot I = 0.18,0.23,0.23 for the no-leak, 
constant-leak and Gaussian leak case, respectively. 

6. Conclusions 

A detailed analysis of NGST multiband photometry applied to the reference case of the study 
of high-redshift galaxies has been carried out with simulations based on galaxy SEDs derived from 
the currently available empirical and model templates and on plausible standard filter-sets. 

In order to correctly identify star forming galaxies in the redshift range 5 < 2; < 9 and early- 
type galaxies above z > 2 and avoid confusion with other SEDs it is mandatory to have photometric 
information in optical bands, besides a standard IR set like FllO F160 K L M. In particular by 
adding the 1^606, /814 and z — Gunn ffiters a good discrimination is obtained above z > 5 for 
star forming galaxies and z > 1 for early-types. The case for an extension of the NGST wavelength 
domain to the optical range is therefore strongly supported by this analysis. 

The effects of leaks in the filter blocking have also been investigated. In spite of rather pes- 
simistic assumptions (a constant leak at a level of 10~^ of the peak transmission over the whole 
spectral range or a leak of Gaussian shape placed at 1.5 times the effective wavelength of the filter 
with an amplitude of 10^^ of the peak transmission and a width of 5% of the central wavelength) 
the effects are not dramatic: the accuracy of the determination of photometric redshifts with a 
standard 1^606 J814 z — Gunn FllO F160 K L M filter set is not significantly deteriorated for a 
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sample limited to Mab = 28. In the range 5.5 < z < 12.5 a rather good accuracy a typical error 
\zsim - Zphotl < 0.2 is achieved. 

As standard specifications for leaks are < 10~^ and as any larger leaks could be measured once 

the filters are produced, it does not look as though leaks constitute a serious problem, at least for 
the particular application described here. The interpretation of fields containing highly reddened 
sources may, however, be more seriously compromised. Extended and continuous spectral coverage 
appears to be the driving factor to maximize the scientific output of the NGST camera. 
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